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ABSTRACT 

The Balloon-borne Large Aperture Submillimeter Telescope (BLAST) operated successfully during a 
250-hour flight over Antarctica in December 2006 (BLAST06). As part of the calibration and pointing 
procedures, the red hypergiant star VY CMa was observed and used as the primary calibrator. Details 
of the overall BLAST06 calibration procedure are discussed. The 1-a uncertainty on the absolute 
calibration is accurate to 9.5, 8.7, and 9.2% at the 250, 350, and 500 pm bands, respectively. The 
errors are highly correlated between bands resulting in much lower error for the derived shape of the 
250-500 pm continuum. The overall pointing error is < 5" rms for the 36, 42, and 60" beams. The 
performance of the optics and pointing systems is discussed. 

Subject headings: balloons — submillimeter — telescopes 



1. INTRODUCTION 

The December 2006 flight of the Balloon-borne Large 
Aperture Submillimeter Telescope (BLAST) incorpo- 
rated a 1.8-m parabolic primary mirror and large- format 
bolometer arrays operating at 250, 350, and 500 pm. A 
complete description o f the BLAST instrument can be 
found in lPascale et~aLl (|2008f ). The BLAST bands sam- 
ple the peak of the spectral energy distribution (SED) for 
cool dust (^10K). Astronomical signals at these wave- 
lengths are difficult or impossible to access from even 
the best ground-based sites. As a result, BLAST has 
the ability to conduct unique Galactic and extragalactic 
submillimeter surveys with sub-arcminute resolution and 
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high sensitivity. BLAST'S primary scientific motivations 
are to study the spatial and redshift distribution and 
evolution of high-redshift star-forming galaxies and to 
probe the earliest stages of star formation within Galac- 
tic molecular clouds. 

BLAST conducted a 250-hour flight, launching from 
McMurdo Station, Antarctica on 2006 December 21, 
and landing on the Antarctic Plateau 2007 January 
2 (BLAST06). BLAST flew at an average altitude 
of 38.6 km with diurnal variations between 37.5 and 
39.6 km. Several extragalactic and Galactic fields 
were mapped, including two large (8.7 deg 2 ) and one 
deep, confusion-limited (0.8 deg 2 ) extragalactic fields and 
two large, overlapping regions (a 50 deg 2 deep region 
and a 200 de g 2 wide region) in the direction of Vela 
(jDevlin et al.ll200l iNeTterfield et al.ll2009h . 

The primary science goals of the BLAST experiment 
demand an absolute flux calibration accuracy of 5-10% in 
all three BLAST pass-bands. In particular, a target, mi- 
correlated 5% uncertainty is driven by the extra-galactic 
science case to enable precise measurements of colors and 
thereby constrain the bolometric l uminosities and star - 
formation rates of distant galaxies (Hughes et al. 2002). 
In order to achieve these goals, a primary calibration 
source for BLAST with the following properties was re- 
quired: (i) point-like and bright enough to be detected in 
each band with a SNR exceeding 20er in a single map; (ii) 
the absolute (correlated) uncertainty in the SED had to 
be less than 10%; and (iii) the uncorrelated components 
of the uncertainty in the SED (i.e. uncertainties in the 
ratios of flux densities in different BLAST bands) could 
be no greater than 5%. In this paper, we report on the 
calibration and performance of BLAST06, concentrating 
on the differences from th e calibration proce dures used 
in B LAST05, discussed in lTruch et all ()2008l ) (hereafter 
IT08I ). Section [2] outlines the basic reduction steps and 
characterization of BLAST06 data. Section [3] discusses 
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the performance of the warm optics in BLAST06. Sec- 
tion |4] outlines the pointing performance of BLAST06. 
Section [5l the bulk of this paper, describes in detail the 
absolute calibration (from detector Volts to Jy on the 
sky) derived from the primary flux calibrator VY CMa 
(|Hoffmeisterl[T93lt iGuthnick fc Schnellerlll939f K 

2. DATA REDUCTION 

The data reductio n for BLAST is dis c ussed in det ail in 
iPascale et all ([2008f ) . IPatanchon et alJ (pOOl ). and[T0l 
Briefly, the data from BLAST consist of a set of 288 
bolometer time streams, in voltage units, sampled at 
100 Hz. Eighteen of these time streams are diagnos- 
tic channels, useful for removing common mode noise, 
the remaining 270 are coupled to the telescope. These 
bolometer data are first cleaned for post-flight analy- 
sis by being de-spiked and then deconvolved to remove 
the effects of the data acquisition system filters from 
the timestreams. The cleaned data are then combined 
with a post-flight pointing solution (|Pascale et al.ll200"8l ) 
to make maps at each wavelength. The map-making pro- 
cess takes advantage of the multiple detectors, as well as 
significant scan cros s-linking, to minimize s triping due to 
instrumental drifts (jPatanchon et al.1l2008| ). 

The bolometers in each array are corrected for relative 
gains, or flat-fielded, so that multi-bolometer maps can 
be generated. The flat-field corrections are determined 
by making individual maps for each bolometer from a 
single point-source calibrator, in this case, VYCMa (see 
§ l5.ip . The bolometers are also corrected for responsivity 
variations over time by using the signals from a calibra- 
tion lamp in the optics box which was pulsed every 15 
minutes throughout the observations. The resulting sig- 
nal is used to correct any time- varying changes in respon- 
sivity per bolometer. Both the time varying changes and 
the variations of beams across each array are small and 
the amo unt o f variation is comparable to those detailed 
in § 2 of [TO! for BLAST05. 

To calculate the flux density from a point source 
we adopt a matched-filtering technique similar to that 
used to extract point sources from s everal recent extra- 
galactic submillim eter surveys Ce.g.. iCoppin et af1l2006t 
iScott et al.|[2006l ) and detailed in § 2 of|T08j. The beam 
profile on the sky, or point spread function (PSF), used 
for calibration and flux extraction is generated by stack- 
ing and averaging several observations of VYCMa af- 
ter removing the ba ckgro und. The PSF is normalized 
as m Equation 1 of pTOl if maps are further filtered, 
the PSF used for analysis of that map must also be fil- 
tered after normalization. See § 15.11 for discussion of the 
background-subtraction technique. 

Beca use the BLAST filt ers have wide (30%) spectral 
bands (jPascale et al.ll2008L § 2.5 and Figure 2), the de- 
rived flux densities depend on the shape of the source 
spectrum within the band, and thus on the temperature 
of the source. A correction must be made to calculate 
m onoc hromatic flux densities from the maps. See § 5.1 
of 1X081 for details of this color-correction. 

3. BLAST06 WARM OPTICS PERFORMANCE 

The BLAST06 warm optics (primary and secondary 
mirrors) performed within their specifications. The 
BLAST06 optical performance was very much improved 
over that of BLAST05 through the use of a new alu- 



minum primary mirror as well as an in-flight focusing 
system ([Pascale et al.ll2008l § 2.4). PSFs for each of the 
BLAST bands are shown in Figure [TJ 

With the pre-flight predicted Noise Equivalent Flux 
Density (NEFD) of 220 mJys 1/2 , we would expect 1- 
a surface brightness fluctuations at the nominal resolu- 
tions of -NEFD/fi, or 11, 4.7, and 2.6 MJy sr" 1 s 1 / 2 at 
250, 350, and 500 pi, respectively where is the area of 
the beam (FWHM 2 ). The measured sensitivities in each 
band were 8.8, 4.8, and 2.7MJysr~ 1 s 1 / 2 , respectively, 
which demonstrates that both the detectors and optical 
efficiencies largely achieved our design goals. The only 
exception is the 250 um beam which is slightly larger than 
the diffraction limit, and has significant side- lobes. For 
each PSF, 76%, 92%, and > 95% of the power is located 
within a Gaussian fit to the PSF, which have FWHM of 
36, 42, and 60" at 250, 350, and 500um, respec t ively. See 
Table [[J See § 3 and Figure 4 of IPascale etaD ([20081 ) for 
further details on the noise and sensitivity of the BLAST 
detectors. 

4. BLAST06 POINTING PERFORMANCE 

Pointing is measured in-flight to an accuracy of ~ 
30" rms using a combination of fine and coarse sen- 
sors, including fiber-optic gyroscopes, optical star cam- 
eras, a 1 _d2ffe£ential GPS, magnetometer, and Sun sen- 
sor (jPascale et al.l [20081 § 7). Post-flight pointing re- 
construction uses only the gyroscopes and day-time star 
cameras. The algorithm is based on a similar multiplica- 
tive extended Kalman filter tec hnique used by WMAP 
(|Pittelkaul feOOlt iMarklevi [20031 ) . modified to allow the 
evaluation of the alig nment parameters o f the star cam- 
eras and gyroscopes (jPascale et al.l [20081 § 11). The off- 
set between the star cameras and the submillimeter tele- 
scope was measured by repeated observations of pointing 
calibrators throughout the flight. We find that the rela- 
tive pointing between the star cameras and submillimeter 
telescope varies as a function of telescope elevation. We 
apply an elevation-dependent correction to pitch and yaw 
with approximate peak-to-peak amplitudes of 260" and 
~ 36", respectively, over the full 25-60° elevation range 
of the telescope. 

Post-flight pointing accuracy is verified by a stacking 
analysis on one of the extragalactic maps. Using the deep 
radio Extended Chandra Deep Field — Sout h Very Large 
Array (E-CDF-S VLA) survey at 1.4 GHz ([Miller et all 
12008! ) we stack patches of the BLAST maps centered at 
the radio source coordinates, summing the flux density 
pixel by pixel (see Figured]). We find that the peak in the 
stacked map is located within 2" from the nominal posi- 
tion of the catalog, more than 15 times smaller than the 
beam size. Moreover, assuming random Gaussian point- 
ing errors, we superimpose the synthetic scaled PSF on 
the stacked map and convolve it with a Gaussian pro- 
file, modeling the broadening of the PSF due to pointing 
jitter. By varying the jitter width, we compute the % 2 
of the convolved PSF over the stacked data. In this way 
we estimate the upper limit in potential random pointing 
errors to be < 5". 

5. ASTRONOMICAL FLUX CALIBRATION 

The primary scientific goals of the BLAST experiment 
demand an absolute flux calibration accuracy of ~10% 
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Fig. 1. — Point Spread Functions (PSFs), provided by observations of VYCMa (S I5.lt . for each of the three wavebands from BLAST06, 
250, 350, and 500 nm f rom left to right. These are generated by stacking several point source maps in telescope coordinates and removing 
the backgro und (§ | 5.2[l. The sm all circles represent the expected diffraction limited FWHM for each of the wavebands, 30, 42, and 60", 
respectively JPascale et alj20 08l. Fitting Gaussians to the PSFs results in FWHMs of 36, 42, and 60" which contain 76%, 92%, and > 95% 
of the power, respectively, indicating the relative power in the sidelobes, which are most visible in the 250 pm PSF. Each PSF has been 
set to zero outside a radius much greater than its FWHM (1.'5, If 5, and 2.'0, respectively) such that > 98% of the power is represented. 
Cutting at larger radii begins to include residual noise at larger scales not removed by the background template. 



TABLE 1 

Calibration and Performance Parameters for BLAST06 



Band 


calib. coeff. 


uncertainty 3 




Pearson correlation matrix b 




PSF FWHM 


[pm] 


[xlO 12 Jy V" 1 ] 


[%) 


250 pm 


350 pm 


500 pm 


[arcsec] 


250 


2.73 


9.5 


1 


0.83 


0.80 


36 


350 


2.86 


8.7 




1 


0.83 


42 


500 


1.16 


9.2 






1 


60 



Note. — Calibration coefficients, calibration uncertainties, and Pearson correlation matrix, showing the relationship between errors in 
different bands, for BLAST06. Note the high correlation means that the measurements involving the ratio of BLAST brightnesses, such 
as spectral indices and temperature, can be reported with much higher accuracy. Also included are the measured FWHM of each beam 
determined by fitting a Gaussian; l-cr errors on the FWHM are of order 1". 

a These values include the 5% uncertainty estimated for the bandpass measurements. Without them, they are 8.1%, 7.1%, and 7.8% at 
250, 350, and 500 pm, respectively. 

b These values include the 5% uncertainty correlated at 50% between bands for the bandpass measurments. Without them, they are 0.977 
between 250 and 350 pm, 0.917 between 250 and 500 pm, and 0.981 between 350 and 500 pm. 



in all three BLAST pass-bands. Achieving this was com- 
plicated by the variability in what part of the sky was 
visible to BLAST due to the unstable projected-latitude 
of the telescope gondola during the flight, the restrictions 
on visibility due to the Sun and Moon avoidance crite- 
ria, the orientation of Sun-shields and other baffling, and 
the elevation r ange ( 25-60°) of the gondola's inner-frame 
(jPascale et all 12003 ). Consequently BLAST had only 
limited access to the calibration sources commonly used 
at submillimeter and far-infrared (FIR) wavelengths. 

Since the Ecliptic plane was not visible during the 
BLAST06 flight, no absolute flux calibration could be de- 
termined from observations of Uranus or Mars, for which 
model SEDs are known to have syste matic uncertain- 
ties < 5% at submil limeter wavelengths ([Griffin fc Ortonl 
[T99llWrightl [2007). The pre-flight strategy for achiev- 
ing a 10% calibration accuracy required the identifica- 
tion of alternate Galactic and extragalactic sources that 
could act as primary and secondary calibrators. The 
requirements included: (i) availability throughout the 
flight; (ii) considered, in some cases, as secondary calibra- 
tors for ground-based submillimeter telescopes and FIR 
satellites; (iii) well-constrained SEDs in the FIR to mm- 
wavelength regime, enabling accurate interpolation of the 



band-averaged flux densities at BLAST wavelengths; (iv) 
bright (^> 1 Jy at 500 pm) and compact sources (with re- 
spect to the BLAST beam-size, i.e. < 20") that reside 
in regions with minimal spatial-structure in the Galactic 
foregrounds or backgrounds, allowing accurate subtrac- 
tion of any extended emission. Given these criteria we 
scheduled regular observations throughout the flight of 
bright embedded protostellar-sources and compact H II 
regions, as well as the star VYCMa. For calibration pur- 
poses, we assume that the SEDs are smooth and neglect 
a ny con tribution from molecular emission lines (See § 4.1 
of fTOg) . 

5.1. VYCMa — The Primary Calibrator for BLAST06 

We chose VYCMa, a red hypergiant star, as our pri- 
mary calibrator. It was the most isolated, point-like 
bright object available during the flight. It is located 
in the Galactic Plane at RA 7 fi 22 m 58^3, Dec -25°46'3"2 
(J2000) and is one of the most intrinsically lumin ous stars 
known ^Humphreys et aT1l2007tlChoi et al.ll2008l ). Unfor- 
tunately, it is not entirely point-like due to the Galactic 
cirrus, which significantly complicates calibration. Fig- 
ure [3] shows the BLAST VYCMa map. An accurate 
PSF is required to use the calibration method as out- 
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Fig. 2. — A cut through the peak-normalized stacked BLAST 
250 pm map at the positions of VLA 1.4 GHz radio sources (dashed 
line) and through the 250 pm PSF (solid line). We see that the 
stack is very well described by the PSF, in both position and width. 
We conclude that our absolute pointing is good to 2" and that 
random pointing errors are < 5" rms. 
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Fig. 3.— BLAST map of VY CMa at 250 pm. This map is in sky 
coordinates and does not have the background removed. 



lined above and in T08. The procedure used to remove 
the background around VYCMa is discussed below, in 



Data are collated to generate the FIR SED used to 
calibrate BLAST06 in a manner similar to Arp 220 from 
BLAST05 (|T0l § 4.1) (see Table H and the SED in 
Figure |4]). These data consist of measurements from 
IRAS at 12, 25, 60, and 100 um (IHelou fc Walker! Il988h . 
from SCUBA on the JCMT at 850 um (T. Jenness, pri- 
vate communication), from SHARC-II on the CSO at 
350pm (D. Dowell, private communication), from Bolo- 
cam on the CSO at 1.1mm (J. Aguirre, private com- 
municat ion), and from UKT14 on the JCMT at 450 and 
800 pm (Sandell 1994|). Careful attention has been placed 
on calculating the full corr elation matrix of the errors, 
as discussed in § 4.1 of T08. Briefly, measurements from 
the same instrument are taken to be 100% correlated (in 
addition to statistically independent photometric uncer- 




1000 100 10 

Wavelength [fim] 

Fig. 4.— Spectral Energy Distribution (SED) of VYCMa, the 
absolute flux calibrator for BLAST06. The best-fit model (heavy 
solid-line) is cons trained using the published data (blue crosses, 
discussed in §[5j] listed in Tabled, excluding BLAST06 measure- 
ments. The grey-lines show the 68% confidence interval, estimated 
from 100 Monte-Carlo simulations, about the best-fit model. Black 
diamonds indicate model predictions for BLAST06 at 250, 350, and 
500 pm. The l-cr uncertainties associated with these predictions are 
8.1%, 7.1%, and 7.8% in each band, respectively, which are highly 
correlated. For reference, the three BLAST passbands are shown 
(normalized to an arbitrary amplitude on the plot). 



tainties), and all measurements are assumed to have an 
additional fully correlated error of 5%, since the bulk of 
these instruments were all ca librated to the same Uranus 
SED (IGriffin fc Ortonl[l99l . 

We find a best-fit model with parameters T = 346 ± 
19 K, /3 = 0.545 ± 0.046, and S Fm = (3.42 ± 0.38) x 
10 17 Wm- 2 (see Figure 0}. This results in BLAST cal- 
ibration uncertainties of 8.1, 7.1 and 7.8% at 250, 350, 
and 500 pm, respectively. These are highly correlated, 
with Pearson correlation coefficients of P250-350 = 0.977, 
P250-500 = 0.917, and P350-500 = 0.981. See Table [TJ 
The large degree of correlation in these values shows that 
almost all of the uncertainty in the SED is in the abso- 
lute value. The relative calibration uncertainty between 
different BLAST bands (i.e. color uncertainty), on the 
other hand, is actually very small, and exceeds the re- 
quirements for the experiment. In fact, care must be 
taken when fitting simple model SEDs to high-SNR mea- 
surements where the calibration uncertainties dominate. 
In such cases, the BLAST data may place such strin- 
gent constraints on SED shapes (e.g. slope and curva- 
ture through the BLAST bands) that single-temperature 
thermal SEDs are not consistent with the observations. 
Section 4.5 of lWiebe et al.l |2009) discusses how to incor- 
porate the bandpass uncertainties for high-SNR sources. 

We include an additional calibration uncertainty re- 
lated to our measurement of the BLAST spectral band- 
passes. Since we have used our measured bandpasses 
to calculate the incident power on our detectors given 
the SED of our primary calibrator, any sources with dif- 
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TABLE 2 

SlJBMILLIMETER FLUX DENSITIES FOR VY CMA 



Wavelength 


Flux Density 


Error 


Instrument 


Reference 


[urn] 


[Jy] 


[Jy] 






12 


9919 


992 


IRAS 


(Helou & Walker 1988) 


25 


6651 


665 


IRAS 


(Helou fc Waiker 1988) 


60 


1453 


115 


IRAS 


(Helou fc Walker 1988) 


100 


331 


33 


IRAS 


(Helou fc Walker 1988) 


350 


15.1 


1.5 


SHARC II 


(D. Dowell, priv. comm.) 


450 


9.7 


4.9 


UKT 14 


(Sandell 1994) 


800 


2.81 


0.4 


UKT 14 


(Sandell 1994) 


850 


2.13 


0.3 


SCUBA 


(T. Jenness, priv. comm.) 


1100 


0.75 


0.08 


BOLOCAM 


(J. Aguirre, priv. comm.) 


250 a 


37.4 


3.7 


BLAST 


This paper 


350 a 


15.0 


1.7 


BLAST 


This paper 


500 a 


6.66 


0.9 


BLAST 


This paper 



Note. — Flux densities used to generate the VYCMa SED shown in Figure [4] as discussed in the text. Flux densities for BLAST and 
from IRAS have been color-corrected. The error column shows photometric and calibration uncertainties added in quadrature for each 
non-BLAST measurement; only calibration uncertainty is shown for BLAST. 
a Extracted from the fit and used to calibrate BLAST06. 



ferent SEDs could have uncertainties in their reported 
band-averaged flux densities arising from uncertainties 
in th e BLAST bandpass measurements (See Eq. 5 of 
IT08D . The dominant source of uncertainty in this mea- 
surement is our kno wledge of the locati on of the band 
edges (See Fig. 2 of iPascale et al.l 12008). which we be- 
lieve to be no greater than about 5%. Assuming this 
maximum uncertainty in the measurement, we have cal- 
culated the effect it would have on sources with different 
SEDs, ranging from "hot" sources (T >40K), to "cool" 
sources (T «15K). Since all of the BLAST filters sam- 
ple the Rayleigh- Jeans regime of VY CMa, the SEDs of 
"hot" sources resemble the primary calibrator and this 
uncertainty has no effect. However, for "cool" sources, 
the SED turns over, particularly in the 250 um band. 
We, therefore, include a 5% uncertainty due to the un- 
certainty in the knowledge of the bandpasses. Since the 
bandpasses were measured simultaneously with the same 
Fourier transform spectrometer, there may be significant 
band-band correlations in the bandpass measurements. 
We assume correlation coefficients of 50% between all 
bands. 

5.2. VYCMa Background Subtraction 

A complication with the use of VYCMa as a calibrator 
is that it is not entirely point-like due to Galactic cirrus. 
We attempt to use the cross-linking (limited to an angle 
between scans of only 9° due to the extreme southerly 
flight of BLAST06) of the region in order to separate the 
beam shape from the diffuse background emission. The 
former is coherent in telescope coordinates, while the lat- 
ter is coherent in sky coordinates. By iterating between 
the two coordinate systems, it is possible to better sep- 
arate the two components, since features in the sky are 
smoothed out to lower significance levels in telescope co- 
ordinates and features in the PSF are smoothed out to 
lower significance in sky coordinates. 

A zeroth-order background subtracted map is gener- 
ated by extrapolating the surrounding emission into the 
region of the point source. A low-order 2-dimensional 
polynomial is fitted to the map with the source masked 
out. 

We start with the true data stream, d, and an initial es- 



timate of the PSF, Pq, taken from the zeroth-order back- 
ground subtracted maps. From there, we scan the esti- 
mate of the point source (PSF), Pj, using a sky simulator 
(an inverse map-maker which takes maps and the BLAST 
pointing solution and generates raw timestreams). The 
iteration index is denoted by i. The PSF maps are in 
telescope coordinates, so scanning must also be done in 
telescope coordinates. The result is a simulated PSF- 
only data stream, p*. We then subtract p* from the 
true data stream, d, producing the background-only data 
stream, bi — d — p* . Then we create a point-source sub- 
tracted background map, B t , from the background-only 
data stream. This map is in sky coordinates. The central 
region of the background-only map, Bi, can be option- 
ally smoothed, but this is only done in the first iteration. 
Next, the estimate of background, Bi, is scanned with 
the sky simulator, producing the simulated data stream, 
b* . Then, the background-only data stream is subtracted 
from the true data stream, producing the PSF-only data 
stream, pi — d — b*. Finally, a new estimate of the PSF, 
Pi+i, is made using the simulated data stream 

We performed 10 iterations of this algorithm in all 
three bands. We examined the RMS in the beam dif- 
ference maps between iterations, finding that it was less 
than 1% of the RMS in the total beam map for the final 
iteration. In addition, the ratio of these quantities was 
shown to drop monotonically, with most of the decrease 
from ^20% to 1% occuring in the first 2-5 iterations, af- 
ter which it flattened out, demonstrating that the answer 
had converged. 

We note that the most complicated beam pattern is at 
250 pm. However, the final answer is least sensitive to the 
details of the backround subtraction at this wavelength 
because the point-source contrast increases. Since the 
beam is smaller, the brightness in a beam of the diffuse 
background is reduced. Furthermore, the diffuse dust in 
the ISM has a temperature of about 20 K (Schlcgel et al. 
I1998L e.g.), whereas VYCMa has a much warmer effec- 
tive SED in the BLAST bands of 200 K, and is therefore 
much brighter at 250 pm. 



6. CONCLUSIONS 
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For the BLAST06 flight, improved optics over the 
BLAST05 flight yielded improved point-source sensitiv- 
ity. BLAST06 calibration uncertainties are approxi- 
mately 10%, as expected. Since we calibrate all BLAST 
bands to one object (VY CMa) the calibration uncertain- 
ties are highly correlated; this high correlation should be 
taken into account when using BLAST photometry mea- 
surements to constrain SEDs. 

Our absolute post-flight reconstructed pointing has 
been shown to be good to 2" with random pointing er- 
rors < 5". Together with the fact that BLAST06 uncer- 
tainties are typically much smaller than those of other 
experimental measurements in the 250-500 pm wave- 
band this should make BLAST data useful as a cross- 
correlation check for future submillimeter experiments. 
Future observation s of VY CMa with t he SPIRE instru- 
ment on Herschel (| Griffin e t al. 2004) would provide a 



useful cross-check of the calibration here, and the higher 
resolution would improve the background subtraction. 
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